Analysis of Localization Phenomena in Weakly Interacting Disordered Lattice Gases 
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Disorder plays a crucial role in many systems particularly in solid state physics. However, the 
disorder in a particular system can usually not be chosen or controlled. We show that the unique 
control available for ultracold atomic gases may be used for the production and observation of 
disordered quantum degenerate gases. A detailed analysis of localization effects for two possible re- 
alizations of a disordered potential is presented. In a theoretical analysis clear localization effects are 
observed when a superlattice is used to provide a quasiperiodic disorder. The effects of localization 
are analyzed by investigating the superfluid fraction and the localization length within the system. 
The theoretical analysis in this paper paves a clear path for the future observation of Anderson-like 
localization in disordered quantum gases. 
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In condensed matter or statistical physics disorder is 
typically neither avoidable nor controllable. It is known 
to lead to spectacular phenomena, such as Anderson lo- 
calization Q, Q. In this case the interference of waves 
scattered on random impurities or defects is responsible 
for localization. In the so called weak localization regime 
coherent backscattering may be regarded as a precursor 
for Anderson localization For a recent treatment of 
light waves scattered by atomic ensembles see e.g. Q. 

Cold atomic gases in magnetic and optical potentials 
offer the unique possibility to introduce well controlled 
disorder to the system. For this purpose several meth- 
ods have been proposed to create disordered, or quasi- 
disordered potentials. These proposals include the use 
of speckle radiation 1^ 0, Q, incommensurable optical 
lattices 0, |MM 0, and the interaction with impu- 
rity atoms [121 . The unique control provided by these 
methods should enable novel studies of disorder induced 
effects, inaccessible to other systems. In addition disor- 
der appears also naturally close to the surface of atom 
chips [l3| and leads to a fragmentation of Bosc-Einstein 
condensates (BEC). 

The first experimental attempts to systematically 
study the effects of disorder on BEC used a speckle po- 
tential. The Florence group 0| observed the apparent 
fragmentation of the condensate into pieces in the pres- 
ence of such a random potential, resulting in character- 
istic stripes in the time-of-flight density profiles. An- 
other experiment, performed simultaneously in Florence 
[lij and in Orsay 0| showed that in the presence of 
the random potential the expansion of the condensate 
is strongly inhibited. A plausible scenario for this sup- 
pression of transport given in [lil IT^ ITsj notes that the 
transport stops when the BEC encounters a random po- 
tential modulation of sufficient height. This explanation 
is purely classical but it is supported by simulations of 



the Gross-Pitaevskii equation (CPE). 

Another disorder configuration was realized in our 
group by combining a one dimensional optical lattice and 
a random potential |l9l | . This potential configuration was 
designed to mimic the situation first envisaged by Ander- 
son However, it was shown that, contrary to first ex- 
pectations, one can not observe Anderson localization in 
this experiment due to interaction effects and the length 
scale of the employed disorder potential. We discuss this 
situation in more detail in the following sections. 

The onset of the Bose-glass phase has recently been 
observed |2Cll | for ultracold atoms in the strongly inter- 
acting regime as suggested by 0- For more details on 
strongly interacting bosonic systems we refer the reader 
to the seminal discussion in [21| . We also note that mix- 
tures of fermions and bosons in random optical potentials 
are of considerable interest . One component in such 
a mixture may in fact serve as the disorder in the system 
fi^l l. Here we concentrate solely on the weakly interact- 
ing case for bosons and investigate the circumstances for 
the appearance of Anderson-like localization. 

The paper is organized as follows. In the first part 
we describe the experimental realization of a weakly in- 
teracting lattice gas with a superimposed slowly varying 
disorder potential and discuss the ground state properties 
of such a system. We observe a classical fragmentation 
of the condensate which is confirmed by a theoretical 
analysis based on the GPE but no signature of Anderson 
localization is found. 

Therefore we theoretically investigate another type of 
disorder which has a shorter correlation length in the sec- 
ond part. Following earher suggestions QiIaEl we con- 
sider the ground state of the condensate in the presence 
of a pseudorandom potential introduced by additional 
optical lattices of different wavelengths. Such a potential 
is called a superlattice We discuss the influence 
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of nonlinear interactions on the occurrence of localization 
phenomena in this system. 

DISORDERED LATTICE POTENTIAL 
Experimental results 

Despite a variety of possible realizations of disorder, 
one method is particularly straightforward to implement 
in current experimental setups. It consists of project- 
ing a disordered optical dipole potential onto the atomic 
sample. Following the theoretical suggestion such re- 
alizations have recently been used to investigate the ef- 
fects of disorder on the ground state Il4l and on the 
dynamics of weakly interacting BECs^4l. Iisl ll^ . 

We perform our experiments with ®^Rb Bose-Einstein 
condensates in an elongated magnetic trap. After laser 
cooling and trapping cold atoms are loaded into a clover- 
leaf magnetic trap with axial and radial frequencies of 
cja: = Stt X 14 Hz and uj± = 27r x 200 Hz, respectively. 
These atoms are evaporatively cooled to quantum de- 
generacy resulting in a final number of condensed atoms 
between 1.5 • 10" and 8 • 10*. 

To carry out experiments in a disordered lattice config- 
uration we use two optical dipole potentials, a ID lattice 
and the disordered potential. Both are derived from a 
Ti:Sa laser operating at a wavelength of A = 825 nm. 
Acousto-optic modulators are used to control the inten- 
sity of each dipole potential and optical fibres are em- 
ployed to deliver the light beams to the experiment. 

The optical lattice is created by retro-reflection of a 
laser beam along the axial direction of the magnetic trap. 
The depth of the optical lattice is typically set to 6.5 Er- 
The recoil energy Er is given by E,. = h^k^/2m, where 
TO denotes the atomic mass and k corresponds to the 
wave number of the optical lattice. The detection system 
is used to image the axial position of the atomic cloud 
as well as the the beam waist. This allows for precise 
positioning of the beam waist with respect to the atomic 
sample. 

The disorder potential is created by illuminating a ran- 
domly structured chrome substrate with a laser beam. 
The partially transmitting substrate is imaged onto the 
atomic sample from a direction orthogonal to the optical 
lattice. Due to the resolution of the imaging system, the 
minimal structure size of the dipole potential is limited to 
a few microns. Figure ^ shows an image of the disorder 
potential at the position of the atomic cloud. We de- 
fine the depth of the disorder potential Va as twice the 
standard deviation of the dipole potential, analogously 
to • Figure El shows the calculated auto correlation 
function 
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FIG. 1: Typical intensity distribution of the disordered opti- 
cal dipole potential. An image of the intensity distribution is 
shown at the top. The intensity variation in a small region 
corresponding to the size of the sample is shown below. 



for a typical realization of the disordered potential. In 
the brackets denote the average over position and 
Vdis{x) represents the disordered potential. Figure [21 
shows that the correlation length of the potential decays 
on a typical length scale of 10 /ito. 
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FIG. 2: Auto correlation function for a typical realization of 
the disordered potential. 

After the production of the BEC our experiments are 
performed as follows. The optical lattice is ramped to its 
final depth within 60 ms. Subsequently the disorder po- 
tential is increased to its final depth within another 60 ms 
followed by a hold time of 20 ms. Finally all potentials 
are switched off and the atomic density distribution is 
measured after 20 ms of free expansion using absorption 
imaging. 

In the expanded atomic density distributions we ob- 
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FIG. 3: Standard deviation a of the observed axial density 
distribution from an inverted parabola after 20 ms of free 
expansion, depending on the disorder strength. The red line 
is a linear fit to the data which serves as a guide to the eye. 



serve significant irregular modulations which depend on 
the strength of the disordered potential. To quantify 
these modulations we fit the central peak of the ex- 
panded lattice gas with an inverted parabolic distribu- 
tion and calculate the standard deviation a of the mea- 
sured density profile from this fit. Figure |21 shows the 
dependence of the standard deviation on the depth of 
the disordered potential. Clearly the deviation from the 
unperturbed parabolic distribution grows with increased 
disorder strength. 

To further investigate the expansion of the disordered 
lattice gas we extract the axial size of the central momen- 
tum peak from the fits of the density with an inverted 
parabolic distribution. Figurc0]shows the resulting sizes 
for four different configurations of the combined potential 
as a function of the atom number. 

The red curve shows the theoretical prediction based 
on the Thomas-Fermi approximation and subsequent 
self-similar expansion for a confinement in the mag- 
netic trap only. The blue curve shows an estimate for the 
expansion of a lattice gas in the absence of disorder. It 
is obtained by using the increased axial size of the cloud 
due to the presence of the optical lattice [i^ as a starting 
point for the self similar expansion. 

The addition of the disorder potential yields a more 
surprising result. Despite its very small potential depth, 
the additional axial confinement due to the disorder po- 
tential leads to a significant increase of the axial size af- 
ter expansion. The theoretical analysis below shows that 
the disordered potential induces strong deviations of the 
density profile from the parabolic envelope. Therefore 
significant deviations from the self similar expansion are 
expected for disordered gases. According to figurcOthcse 
deviations produce irregular density modulations and can 
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FIG. 4: Size of the central peak after 20 ms of free expansion 
versus the number of atoms. The clouds are released from the 
following potentials: magnetic trap (red o), magnetic trap and 
disorder (black □), magnetic trap and optical lattice (blue <), 
magnetic trap and disorder and lattice (green l>). The lines 
correspond to a theoretical prediction (see text). The lattice 
depth is 6.5 Er and the disorder has a depth of 0.1 Er. 



lead to pronounced changes in the widths of the expanded 
clouds. The explicit form of the expanded density distri- 
bution however depends strongly on the exact realization 
of the disordered potential. We have used a 3D numer- 
ical simulation to check the expansion of the disordered 
gas in the absence of the optical lattice. Depending on 
the disorder potential used, the simulation confirms the 
observed behavior qualitatively. 



Theoretical Analysis 

In the system considered here, the disorder potential is 
imposed along one axis of the cloud. Hence a ID Gross- 
Pitaevskii equation is used to describe the basic prop- 
erties of the system and to analyze parameter regimes 
where localization phenomena can be observed. The 
GPE provides an appropriate description of the system 
as long as depletion effects are small. We have estimated 
the number of atoms depleted from the condensate wave- 
function within the Bogoliubov theory for the optical lat- 
tice potential used in the experiment. For a total particle 
number of 10^ atoms the fraction of depleted atoms is less 
than 1%. 

The ID GPE used in our simulations is given by 



idt(t) - 
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Vocos^{kx) + VdUx)+gW 



(2) 

where we have adopted harmonic oscillator units, i.e. 
Tiujx, \/'hlmu)x and l/w^; as energy, length and time 
units. The depth of the optical lattice is given by Vb 
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FIG. 5: Top panel: Ground state solution of the GPE for 
a condensate in the combined potential of the magnetic trap 
{ujx = 27r X 14 Hz), optical lattice and disorder potential. Bot- 
tom panel: Central part of the ground state solution of the 
GPE (solid line) and the corresponding solution within the 
effective mass approach (dashed line) . The depth of the opti- 
cal lattice is 6.5 Er while the depth of the disorder potential 
0.7 Er. 



and Vdis(a;) denotes the disordered potential. In the 3D 
case the coupling constant is given by g^D = Anfi^a/m, 
where a denotes the s-wave scattering length. For the ID 
simulations the coupling constant g is chosen such that 
the Thomas-Fermi radius of the ID gas equals the axial 
Thomas-Fermi radius in the 3D trap. 

To simulate the experimental situation the wavelength 
of the lattice is set to A = 825 nm. Since the harmonic 
trap and the disorder potential change on a length scale 
much greater than the lattice spacing and the condensate 
healing length, / = l/\/87ma (where n is the condensate 
density) we eliminate the lattice potential in the GPE by 
applying the so called effective mass analysis [2^ . 

Within this analysis the GPE (0) is replaced by an 
equation where the optical lattice potential is absent but 
the mass of a particle and the coupling constant are cor- 
rected. We assume that the ground state solution of the 
GPE has the form 



2m* 



Y 



+ Vais{x)+g*\f{xf 



(4) 

where m* and g* are the effective mass and the effec- 
tive coupling constant, respectively. For an optical lattice 
depth of 6.5 Er the effective parameters are 



m* = 2.56 m, 
g* = \mg. 



(5) 



For the total number of = 10^ atoms (that implies 
g = 1800) the values of the effective parameters suggest 
that one can employ the Thomas-Fermi approximation 
and neglect the kinetic energy term. Then the solution 
of Eq. ^ is 



2 ^ ^l* -x'^ji- VAi,{x) 
9* 



(6) 



where ^* is determined from the normalization condi- 
tion 



\f{x)\'dx = l. 



(7) 



In Fig. we show a comparison of the ground state solu- 
tion of the full GPE 121) and the solution of the effective 
mass approach obtained within the Thomas-Fermi ap- 
proximation. The squared overlap of these solutions is 
greater than 0.99 and hence they are practically identi- 
cal. 

This analysis shows that the effect of a slowly vary- 
ing disorder potential applied to Bose-Einstein conden- 
sates in the lattice potential can be described within the 
Thomas-Fermi approximation. The condensate density 
is modulated by the slowly varying disorder and conse- 
quently no Anderson localization is present. 

These results thus bear similarity to the experiments 
performed in the absence of an optical lattice [lj| where 
a fragmentation of the BEC is induced by the disorder 
potential. 



SUPERLATTICE POTENTIAL 



0o(x) = VMf{x)m{x). (3) 

Here uo{x) denotes the Bloch function corresponding 
to the ground state of the Schrodinger equation with the 
optical lattice potential, f{x) is an envelope function and 
is a normalization factor. Substituting the ansatz (O 
into Eq. (j2J) leads to the equation 



The above analysis suggests that it is necessary to in- 
troduce a disorder that changes on a length scale smaller 
than the healing length to enter a regime where localiza- 
tion effects can be observed. To overcome the experimen- 
tal difficulties of imposing such a truly random potential 
the use of pseudorandom potentials has been suggested 

m 

These pseudorandom potentials can be formed by one 
or more additional optical lattices creating a so called 
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FIG. 6: Ground state density of a non-interacting BEC in the 
combined potential of a magnetic trap {ujx = 27r x 4 Hz) and 
an optical lattice with a wavelength of A = 825 nm and a 
depth of 6.5 Er- Note the logarithmic vertical scale. 
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FIG. 7: Ground state of a non-interacting BEC in the po- 
tential of Fig|S|with two additional weak optical lattices at 
A = 960 nm and A = 1060 nm. with depths of 0.2Er. 



superlattice. Since cold atomic gases present a finite sized 
system, a suitably chosen pseudorandom potential can 
provide the desired disorder [27j . 

Consider first a situation without interactions as dis- 
cussed previously in Figure shows the ground state 
density for a condensate in the combined potential of the 
magnetic trap and the optical lattice. As expected the 
modulation of the density due to the lattice is visible. 
The addition of two very weak lattices at wavelengths of 
A = 960nm and A = 1060nm changes this situation dras- 
tically as shown in Fig. [7| The exponential tails of the 
ensemble density are a clear manifestation of Anderson- 
like localization. 

However, it is well known that localization phenomena 
are strongly influenced by the presence of interactions 
j28l |. Therefore a simulation including these interactions 
is necessary to predict parameter regimes for the obser- 
vation of localization phenomena. 




FIG. 8: Top panel: Ground state solution of the GPE for a 
condensate in the combined potential of the magnetic trap, 
optical lattice and disorder potential. Bottom panel: Disorder 
potential Vdis{x) (dashed line) and disorder potential plus the 
g\(f>o{x)f term (solid line). The depth of the optical lattice 
is 6.5 Er while the depth of the disorder potential is 0.4 Er- 
The coupling constant is g = 1800. 



Screening due to interaction 

Figure |S1 shows the ground state solution of the GPE 
for a disordered potential that changes on the length scale 
of the optical lattice. While the effective mass approach 
can no longer be applied in this case, the localization 
signatures still do not emerge. This is due to the fact that 
the interaction between atoms delocalize the condensate. 

The effect can be visualized by considering an effective 
potential. That is, once the solution 0o of the station- 
ary GPE is known, the equation can be considered as a 
Schrodinger equation 



1 d^M^) 



2 dx^ 
with the effective potential 



+ Vcs{x)(l)o{x) = fJ-4>aix), 



(8) 
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x 

Y 



Vq cos2(fcx) + Vdi,{x) + g\(j>o{x)?- (9) 



As shown in Fig. [3 the situation without the nonlin- 
ear term g\(j)Q{x)\'^ leads to the observation of Anderson- 
like localization. However, when this term is present the 
effective potential loses its disordered nature due to a 
screening effect. Figure |H1 shows the potentials Vdis(a:) 
and Vdis(a;) + g\(f)Q{x)\'^ . It is apparent that in the ef- 
fective potential the disorder is smoothed and screened. 
For typical experimental parameters the term g|(/)o(x)p 
dominates over Vdis {x) and consequently the randomness 
necessary for localization is lost. 
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Hence a fine scale disorder alone is not sufficient to in- 
duce a non-trivial localization in the system. Additional 
control of the interaction is necessary. This effect is due 
to the accumulation of atoms in the wells of the random 
potential. In these regions of high density the nonlin- 
ear term in the GPE effectively smoothes the potential 
modulations . A detailed discussion of this screening 
effect was recently given in jsoj . 



Analysis of Localization 

We have investigated the transition from a fully de- 
localized to a localized ground state of a bosonic lattice 
gas in two ways. First a range of interaction strengths 
g\(j>o{x)\'^ were investigated numerically for a fixed depth 
of the additional lattices. In current experiments a varia- 
tion of the interaction strength can be realized by varying 
the density |(/)o(a;)p of the sample or by using a Feshbach 
resonance to change the coupling constant g. Our main 
aim was to investigate if suitable parameters for the ex- 
perimental observation of localization can be obtained. 
Secondly we have investigated the onset of localization 
as the disorder strength is increased in a weakly interact- 
ing BEG. 

Our calculations were performed in an experimentally 
accessible regime with a trap frequency of 27r x 4 Hz and 
a pseudorandom potential equivalent to the one used for 
Fig. 13 For (7 = one obtains Anderson-like localization 
of the ground state wavefunction which is characterized 
by an exponential localization as shown in Fig[71 

Figure ini shows ground states within this potential for 
three values of the interaction parameter. As g is in- 
creased the number of localization centers grows and for 
large values of g they overlap considerably. This behav- 
ior suggests that the condensate wavefunction becomes 
a combination of these localized states due to nonlin- 
ear interactions. When g is of the order of 500 one can 
no longer distinguish individual localized states and the 
clear signature of non-trivial localization vanishes. 

The case of g = 256 shown in Fig. ^ is equivalent to 
a 3D experimental realization with trap frequencies of 
cja; = 27r X 4 Hz and wj, = 27r x 40 Hz and N = 10* atoms. 
The ground state simulation shows that characteristic 
features of Anderson-like localization arc present while 
the experimental parameters are within reach. 

In a second part of our analysis we have investigated 
the localization length as a function of the disorder depth. 
The localization length I was obtained by fitting the indi- 
vidual localization sites within the ground state density 
by \4'o{x)\'^ oc exp(-|a; - xo\/l). 

FigurelTUlshows these localization lengths for depths of 
the superlattice potential up to 0.2 Er- The simulation 
was performed for interaction strengths of g = and 
(7 = 8 in a trap with w^; = 27r x 4 Hz trapping frequency. 
Each point in the figure is the result of a fit to peaks in the 
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FIG. 9: Ground states of the GPE for a condensate in the 
combined potential of the magnetic trap, optical lattice and 
pseudorandom potential. The depth of the optical lattice is 
6.5 Er while the depths of the additional lattices forming the 
pseudorandom potential are 0.2 Er- The coupling constants g 
for the panels are: 0.5 (top), 8 (middle), 256 (bottom). Note 
the different scales in the panels. 
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FIG. 10: Localization length as a function of the depth of 
the additional optical lattice potentials for two interaction 
strengths g — (circles) and g = 8 (full dots). The depth of 
the main optical lattice is 6.5 Er. 



wavefunctions depicted in Fig.lHl Within the scope of our 
simulation the amplitude of the pseudorandom potential 
has to be sufficiently big to ensure that the localization 
length is smaller than the size of the system. 

This analysis clearly shows two distinct features of lo- 
calization. First, the localization length strongly depends 
on the depth of the superlattice. Even a very small 
added pseudorandom potentials leads to localization in 
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FIG. 11: Superfluid fraction as a function of the coupling 
constant g for a fixed superlattice depth of 0.2 Er- Full dots 
are obtained for a box potential corresponding to a trap fre- 
quency of Ljj; = 27r X 14 Hz and open circles toujx = 27r x 4 Hz. 
Full squares show the case without a superlattice potential for 
iOx = 2% X 14 Hz. 



small localization sites. This confirms the expected non- 
perturbative character of localization. Secondly the anal- 
ysis shows that the interaction strength only has a small 
effect on the localization length in the strongly localized 
regime as shown in Fig. 1101 This can also be inferred 
from Fig. ini taking into accomit the different axis scales. 



Analysis of Superfluidity 

To further investigate the effect of interactions we have 
analyzed the superfluid fraction of the sample as a func- 
tion of the coupling constant g in a potential box. 

The superfluid fraction is obtained by calculating the 
response of the condensed sample to twisted boundary 
conditions |33|. Within this model the superfluid frac- 
tion acquires additional kinetic energy due to a spatially 
varying phase. A comparison of the energy with and 
without phase twist yields the superfluid fraction which 
is defined as fs = 2[Eo{v) - Eo{Q)]/Nv^. Here Eo{v) is 
the ground state energy when a velocity field v is im- 
posed on the system (i.e. we compute the ground state 
solution in the form (p^ix) cx.Y){ivx) where (j^oi^) fulfills 
periodic boundary conditions) [sJl. In our calculations 
the size of the potential box was chosen to match the 
size of the atomic cloud in the harmonic potential. It is 
important to note that this method only represents one 
possible definition of superfluidity |^ ■ Within the 
GPE framework, it does not include higher excitation 
modes or atoms depleted from the condensate fioj . 

FigureE]shows the superfluid fraction as a function of 
the coupling constant. In the presence of an optical lat- 
tice some loss of the superfluid fraction is observed, but 
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FIG. 12: Ground states of the GPE (previously shown in 
Fig.^lon a logarithmic scale) for a condensate in the combined 
potential of the magnetic trap, optical lattice and pseudoran- 
dom potential. The coupling constant g for the panels are: 
0.5 (top), 8 (middle), 256 (bottom). 



even for small interaction strengths a considerable super- 
fluid fraction remains. This behavior drastically changes 
when the pseudorandom disorder created by the two ad- 
ditional optical lattices at 960 nm and 1060 nm is added. 
At low values of the coupling constant g the superfluid 
fraction strongly decreases, indicating the onset of local- 
ization. However, our analysis shows that the superfluid 
fraction remains large for the coupling constants g in typ- 
ical experimental realizations, indicating the absence of 
Anderson-like localization. 



Detection of localization 

The observation of features of Anderson-like localiza- 
tion may pose considerable experimental difficulties. Fig- 
ure El shows the density distributions on a linear scale. 
In an experimental measurement the exponential nature 
of the density variation within the localization sites will 
probably not be visible. Most likely, limitations due to 
the imaging optics will inhibit the observation of individ- 
ual localization sites for experimentally accessible densi- 
ties. 

However, the onset of localization leads to a consid- 
erable change in the ground state density when a small 
pseudorandom potentials is added. If the localization 
effect is indeed non-perturbative, it may be possible to 
detect this change of the density even for small added 
pseudorandom potentials. 

A second avenue for the detection of localization is a 
time-of-flight measurement of the velocity distribution. 
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FIG. 13: Atomic density after 20 ms of free expansion for 
a condensate prepared in the states shown in Fig. 1121 (left 
column) and without disorder potential (right column). Os- 
cillator units corresponding to a trap frequency of 27r x 4 Hz 
are used. 



Within the experimental part a realization of a disor- 
dered lattice gas is described in detail. An analysis of 
the time-of-flight signal allows for a conclusive evalua- 
tion of possible localization phenomena in the sample. 
In accordance with the theoretical findings however, no 
localization effects are observed. 

The theoretical work within this paper paves a clear 
path for the future observation of Anderson-like localiza- 
tion in cold atomic samples. 
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The simulated density distribution after 20 ms of free 
expansion is shown in Fig. 1131 Despite a clear difference 
in the ground state wavefunction, the width of the enve- 
lope of the zero-momentum peak is strikingly similar for 
all realizations. The results represent a clear distinction 
from the classical case, where the width of the expanded 
cloud depends strongly on the interaction parameter g. 
We conclude that the width of the zero-momentum peak 
mainly depends on the localization length /, which does 
not vary significantly as a function of g in the simulations 
presented here. 

These two options show a path towards the observation 
of Anderson-like localization in the regime discussed. 



CONCLUSION 

We have presented a detailed theoretical analysis of 
two configurations for the production and observation 
of disordered quantum degenerate gases. This analy- 
sis showed that localization effects are not expected for 
the configurations involving a slowly varying disordered 
dipole potential. However, clear localization effects and a 
reduction of the superfluid fraction were observed when 
a supcrlattice was used to provide the disorder. These 
effects can be suppressed due to screening by nonlinear 
interactions within the sample. The dependence of the lo- 
calization features on these interactions and on the depth 
of the superlattice potential were analyzed in detail. It 
was shown that an analysis of the time-of-flight signal will 
allow for a conclusive evaluation of possible localization 
phenomena in the sample. 
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